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ABSTRACT 

We have studied the maser emission from a thin, planar, gaseous ring in Keplerian rotation around 
a central mass observed edge-on. The absorption coefficient within the ring is assumed to follow a 
power law dependence with the distance from the central mass as, k = Kgr"*. We have calculated 
position-velocity diagrams for the most intense maser features, for different values of the exponent q. 
We have found that, depending on the value of q, these diagrams can be qualitatively different. The 
most intense maser emission at a given velocity can either come mainly from regions close to the inner 
or outer edges of the amplifying ring or from the line perpendicular to the line of sight and passing 
through the central mass (as is commonly assumed). Particularly, when q > 1 the position- velocity 
diagram is qualitatively similar to the one observed for the water maser emission in the nucleus of the 
galaxy NGC 4258. In the context of this simple model, we conclude that in this object the absorption 
coefficient depends on the radius of the amplifying ring as a decreasing function, in order to have 
significant emission coming from the inner edge of the ring. 
Subject headings: galaxies: individual (NGC 4258) — galaxies: nuclei — masers 



1. INTRODUCTION 

The a priori probability of seeing a thin disk nearly 
edge-on is very small. It is given by p ~ 0.125 (/i/i?)^, 
where h is the thickness of the disk and R is its radius. 
Typically h/R ~ 0.01 and thus p ~ 1.25 x 10"^ Sur- 
prisingly however, the maser emission observed in sev- 
eral cosmic sources has been successfully modeled as 
coming from a ring or truncated disk in Keplerian ro- 
tation (around a massive object) seen edge-on. For in- 
stanc e: circumstellar disks in star-for ming regions as in 
S255 (|Cesaroni |[l99l and MWC 349 (iPonomarev et al. I 
Il994[ ). and also circumnuclear disks around black holes 
of galactic nuclei as in NGC 4258 (jWatson fc Wallin I 
119941: iMivoshi et al. |[l995l) . In general, the maser emis- 
sion from a Keplerian d isk observed edge-on produ ces a 
triple-peaked spectrum fElmegreen fc Morris lll9"79h : but 
IPonomarev et al71 (| 19941 showed that, there is a transi- 
tion from triple- to double-peaked spectra as the width 
of the amplifying ring decreases. 

NGC 4258 is a Seyfe rt 2 /LINER located at a dis- 
tance of 7.2 ± 0.3 Mpc (jHerrnstein et al. I [19991 ). The 
water maser emi ssion (22 GHz) toward this galaxy was 
first detected bv IClaussen et al.~l ()1984l ). Shortly after- 
wards it was shown that the water masers are confined 
in a very small region (^1 .3 pc) at the ce nter of NGC 
4258 ("Cla ussen fc Lo1ll986[ ). Subsequentlv. iNakai et all 
11993) discovered water maser emission with velocity off- 



Electronic address: 
Electronic address: 



Luscanga@a£trosmoainam.mx| 



ragaignucleares . unam . mx| 



sets ±1000 km s~^ from the already known emission at 
the galactic systemic velocity of ~472 km s^^. They sug- 
gested that the high-velocity emission could arise from 
masers orbiting a massive central black hole, or ejected 
in a bipol ar outflow. Using the Very Long Baseline Array 
(VLBA). lMiyoshi et al. I (|1995l ) simultaneously observed 
the systemic and high-velocity water maser emission in 
NGC 4258, finding that the spatial distribution and line- 
of-sight velocities of the water masers trace a thin molec- 
ular ring in Keplerian rotation around a massive black 
hole of 3.6x10^ Mq seen nearly edge-on. The position- 
velocity (PV) diagram for the maser emission shows dis- 
tinct Keplerian orbits (with deviations < 1%) defined by 
the high-velocity maser emission that arises on the ring 
diameter perpendicular to the line of sight, as well as a 
line traced by the systemic maser emission that arises 
from material on the inner edge of the amplifying ring, 
this linear dependence is a consequence of the change in 
the line-of-sight projection of the rotation velocity. 

By monitoring both systematic and high-velocity wa- 
ter maser emission of NGC 4258 over periods of sev- 
eral years with different radio telescopes, a signif- 
icant centripetal acceleration was observed only for 
the maser features near the galactic systemic velocity. 
The systemic m aser features drift at a mean rate of 
~9 km s~^vr-^ (iHa schick e t al. 1 119941: [Greenhill et al. I 
119951 INakai et aLT l995: B ragg et al. 1 boOOf ) while the 
high-veloci ty ma. s er fea tures drift by <1 km s~"'^yr~^ 
dGrecnh iU et al. I [19951). In a recent spectroscopic 
study, iBragg et al. I ( 2000( ) detected accelerations for 
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the high- velocity features in the range of —0.77 to 
0.38 km yr~^. These measurements indicate that 
the systemic water masers he within a relatively nar- 
row range of radii, on the near side of the ring at the 
proximity of its inner edge, while the high-velocity wa- 
ter masers are located near the ring diameter (bet ween 
-13?6 and 9?3 of the mid-hne. iB^agg et al. II2000D . In 
addition, the deviation of the high-velocity masers from 
a straight line passing through the systemic masers in 
the plane of the sky suggests tha,t the rotating dis k is 
slightly warpe d (iHerrnstein et aTpil996l [l999l. 120051). 

Previously, iWatson fc Wallini (|1994[ ) demonstrated 
that the maser emission from a rapidly rotating, thin 
Keplerian ring viewed edge-on can reproduce the general 
features of the observed 22 GHz radiation from the nu- 
cleus of NGC 4258, including the high- velocity satellites. 
However, it is important to point out that their assump- 
tion of a uniform absorption coefficient within the ampli- 
fying ring results in a PV diagram for the most intense 
masers that is qualitatively different from the observed 
one. While their model predicts that the maser emission 
at velocities around the systemic velocity of the galaxy 
comes mainly from the outer edge of the ring, the ob- 
servations indicate that this emission is actually coming 
from the inner parts of the truncated disk. 

In this paper we show that this discrepancy can be 
resolved if the absorption coefficient decreases with dis- 
tance from the central mass. The model is presented in 
§2. The main results are described in §3. Finally, the 
conclusions are discussed in §4. 

2. MODEL 

We study the maser emission that arises from a thin, 
planar, gaseous ring in Keplerian rotation around a mas- 
sive central object when it is observed edge-on. The mas- 
ing gas is located between Rq and R, the inner and outer 
radii of the ring, respectively. For simplicity, we assume 
that the disk is transparent to the maser radiation at 
radii smaller than Rq and greater than R, although the 
inner region is probably thermalized due to the higher 
gas density, and actually it would absorb a significant 
fraction of the maser radiation produced in the far side 
of the ring (see §4). The absorption coefficient is assumed 
to follow a power law function of the distance from the 
central mass within the amplifying ring as, k = Kgr"^. 
The distances are measured in units of R, and the veloc- 
ities are measured in units of Vout, the rotation velocity 
at the outer edge of the ring (see Figure 1). 

For the case of an unsaturated maser and neglecting 
the spontaneous emission, the intensity of the maser ra- 
diation from a line of sight with impact parameter y at 
a velocity Vr is 



and 



I{vr,y) ^ Ioe^<^^-^\ 



(1) 



where the optical depth or gain along the line of sight is 
given by 



-{V - Vr^ 



T{vr,v)^2nQ / (a:2+2/2)-?/2 ^xp 
with 

^/rl-y"^ for < |y| < ro , 
for ro < |y| < 1 , 



(2) 



The line-of-sight velocity component of the gas at the 
position (x, y) can be expressed as v = y/{x'^ + y"^)^^^. 
Here Iq and Avd are the background intensity and the 
Doppler width, respectively, which are supposed to be 
uniform inside the amplifying ring. The Doppler width 
Avd is related with the FWHM of the velocity distribu- 
tion of the emitting particles as Avd = FWHM/V4 ln2. 

We have numerically solved equations (1) and (2), and 
we have also calculated the y-positions (impact parame- 
ters) of maximum maser intensity for each specific value 
of the velocity Vj.. When we have found two local max- 
ima, we have kept both. With this information, we have 
constructed the PV diagrams using the positions of the 
observer's line of sight with maximum emission at each 
velocity. This wa y to construct t he PV diagrams was 
previously used bv lUscanga et al. 1 (|2005). 

We show the results using the following values for the 
model parameters which seem to be appropriate for mod- 
eling maser emission in the galaxy NGC 4258. The back- 
ground intensity is Iq = 1.3 x 10~^ Jy beam~^, corre- 
sponding to a r adio continuum source w ith a temper- 
ature of 10^ K (jWatson fc Wallin I [l99l . The dimen- 
sionless inner radius rp = 0.51, using the estimated val- 
ues for the inner and outer radii of 4.1 and 8.0 mas re- 
spectively, given by Miyoshi et al. (1995). The Doppler 
width Avd — 0.007wout which combined with an oute r 
rotation velocity of 770 km s~^ ((Mivoshi et al. |[T995f) . 
gives a Doppler width :^ 5 km s~ ^, similar to the value 
used bv IWatson fc Wallin I (jl994l ). We have used some 
representative values of the exponent q, specifically q = 
0, 1/2, 15/8, 5 for Models I, II, HI and IV, respec- 
tively. In Model I, we study the simplest situation of a 
uniform absorption coefficient. In Model III, we choose 
q = 15/8, that corresponds to the densit y dependence 
with t he radius of an accretion disk, i.e., iFrank et al. I 
(jT992). Finally, in Models II and IV, we explore two 
other different values of the exponent q in order to study 
how it changes the results. In all the models, the value 
for the absorption coefficient kq is mainly determined 
by the requirement that the intensity at the peak of the 
central component (13 Jy beam~^) is compatible with 
the observational data when the background intensity is 
1.3x10"^ Jy beam~^. Other values of Iq, tq, Avd, and 
Kq give qualitatively similar results. 

We present the results in the next section; but let us 
first discuss briefly some important concepts in order to 
understand these results. In general, the observed emis- 
sion at a given velocity coming from a specific position in 
a nebula has contributions of the whole material along 
the line of sight. However, when the velocity gradient 
along the line of sight is greater than the dispersion ve- 
locity (thermal or turbulent) of the emitting material, 
the main contribution to the emission is actually coming 
from a narrow region around the point with a line-of- 
sight velocity equal to the observation velocity. The esti- 
mated width of the region is 21, where I is the correlation 
distance defined as 

\av/ax\ 

here dv/dx is the line-of-sight velocity gradient. In this 
approximation, known as Sobolev's approximation or the 
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approximation of high velocity gradient, the observed in- 
tensity is given by the following expression 

I{vr) = /oK)e-^('''-) + S{vr){l - e-^^"-)) , (4) 

where /q is the background intensity, S is the source func- 
tion and T is the optical depth given by 

T{Vr) = k{21) , (5) 

where k is the absorption coefficient. 

For the case of maser emission, the value of k is in- 
trinsically negative and r is also negative, therefore the 
factor e"'^*^"'^^ becomes an amplification factor. Because 
of this reason, the relative contribution at a given veloc- 
ity of the correlation region is even more important with 
respect to the remainder of the emitting material than 
in the case of non-maser emission. Consequently, the ap- 
proximation given by equations (4) and (5) is suitable 
for maser emission. 

As shown in the next section, for a gaseous ring of in- 
ner radius Rq and outer radius R in Keplerian rotation 
and seen edge-on, the emission either comes preferen- 
tially from the inner or outer edges of the ring or from 
the line perpendicular to the line of sight and passing 
through the ring center. In the first two cases, it is easy 
to show that the expected PV diagram will be a straight 
line. When the emission comes from the outer edge, the 
slope of the straight line is equal to one (measuring the 
distances in units of the outer radius of the ring and the 
velocities in units of the rotation velocity at that point), 
whereas if the emission comes from the inner edge, the 

slope of the straight line is equal to l/r^^"^. On the other 
hand, when the emission arises from the line perpendic- 
ular to the line of sight, the PV diagram will be a curve 
with the form l/y^/"^, where y is the impact parameter 
of the observation (see Figure 2) . 

3. RESULTS 

The PV diagrams for the maser emission peak are 
point-symmetric, consequently we only discuss positive 
velocities from now on (see Figures 3 and 4). 

• Model I (g = 0) - With this value of the expo- 
nent q, we are considering the simplest situation, 
a rmiform or constant absorption coefficient. The 
strongest maser emission cither comes mostly from 
the outer edge of the ring at velocities lower than 
1, or from the mid- line of the ring perpendicular 
to the line of sight and passing through the cen- 
tral mass at greater velocities. The filled squares, 
circles, and triangles mark the regions of strongest 
maser emission at each velocity. 

• Model II {q = 1/2) The results are qualitatively 

similar to those of Model I. 

• Model III {q = 15/8) - This value of exponent q 
corresponds to the density dependence with radius 
of an accretion disk {p cx r^^^/*). The strongest 
maser emission cither comes mainly from the inner 
edge of the ring at low velocities (velocities near 
the systemic velocity), or from the outer edge at 
velocities close to 1. On the other hand, at ve- 
locities greater than 1, the most intense emission 
comes predominantly from the mid-line of the ring 
perpendicular to the line of sight. 



• Model IV {q = b) - The strongest maser emission 
comes mainly from the inner edge of the ring at 
velocities lower than 1. At greater velocities, the 
most intense emission can either come mainly from 
the inner edge or from the mid-line of the ring per- 
pendicular to the line of sight. 

In summary, from the results of Models I -IV (see Fig- 
ure 3), we found that the most intense maser emission 
can be around the inner or outer edges of the ring, or 
the mid-line of the ring perpendicular to the line of sight 
depending on the velocity and also on the value of q. In 
fact, the PV diagrams are qualitatively different when 
q < \ or q > \. In the first case, for q < 1 (including 
the simplest situation with a uniform absorption coeffi- 
cient, g = 0) and < 1, the PV diagram corresponds 
to a straight line with slope 1; for > 1, the diagram 
corresponds to a Keplerian curve. In the second case, 
for q > 1 and Wj. < 1, the PV diagram corresponds to a 
straight line with a slope that depends on the inner ra- 
dius of the ring. At velocities close to 1 the slope changes 
to 1; for > 1, the diagram corresponds to a Keplerian 
curve and also a straight line with a slope that depends 
on the inner radius under circumstances such as in Model 
IV. 

It is also important to realize that when g > 1, the 
optical depth or gain presents two local maxima within 
a certain velocity range. Either local maxima may be a 
global maximum. For Vr < 1, the local maximum can be 
either at the inner and/or outer edges of the ring, while 
for Vr > I, they are located at the inner edge and/or 
mid-line of the ring (see Figure 5). 

As shown in the top panels of Figure 5 (Model III, 
q = 15/8), the relative difference between the two lo- 
cal maxima is not very significant. However, when the 
value of q is higher (like in Model IV, q = 5 shown in 
the bottom panels) the relative difference becomes more 
important. 

In order to estimate the velocity Vc at which the global 
maximum of the optical depth changes its locus, we have 
calculated analytical approximations for the largest value 
of the optical depth or gain that corresponds to the max- 
imum intensity at the inner and outer edges of the ring, 
and also at the mid-line of the ring perpendicular to the 
line of sight. The detailed calculations are presented in 
the Appendix. The following equations give the local 
maximum depth as a function of the velocity in each 
neighborhood 



1-9 
'^0 



4 Tn 
T{Vr) O:^ -VtTKo^Vd 7= ^ 

o Vr^/l-rov:^ 
4 ^ 

T{Vr) ~ -^/tTKqAvd 



inner edge , (6) 

7== outer edge , (7) 

T{vr) — -^Ko\J Au/jiJr'"^ mid-line. (8) 
V 3 

The velocity is estimated by combining equations 

(6) and (7), or equations (6) and (8) according to the 
value of Vc (when < 1 or t^c > 1, respectively). The 
results are 



- I 



2(l-g) 
'^0 



1/2 



for Vc<l, (9) 



4 



USCANGA, CANTO, & RAGA 



^AvDrl^'^"'' - vt"-"" + rov^'-' = for Vc > 1, (10) 

which are presented in Figure 6 for some representative 
values of the exponent q. 

The bottom plot of Figure 6 shows Vc as function of 
the inner radius ro for different values of q, from equation 
(9a). For q < 1, there is no solution to equation (9a). 
When q = 1, Vc = for any value of rg. That is, the 
optical depth has a maximum and its locus is around 
the outer edge of the ring, and Vc is meaningless as we 
have defined it. When q > I, the optical depth presents 
two local maxima, and Vc is different from zero and its 
value depends on rg. This velocity corresponds to the 
value at which the locus of the global maximum changes 
from the inner edge to the outer edge of the ring. As a 
consequence, there is a slope change in the PV diagrams 
at velocities lower than the rotation velocity at the outer 
edge of the ring. For instance, when q = 15/8 and tq = 
0.51, the slope change occurs at Vc — 0.906. In other 
words, the locus of the global maximum of the optical 
depth changes from the inner to the outer edge of the 
ring at this velocity Vc- 

The top plot of Figure 6 also shows Vc as function of 
the inner radius tq using specific values of q and Avd in 
equation (9b); in this case, 5 and O.OOTiiout, respectively. 
As an example, when Avjj — 0.007wout, g = 5 and rg = 
0.51, then Vc — 1.077. Stated differently, at that velocity 
Vc, the largest value of both the optical depth and the 
intensity changes its locus from the inner edge to the 
mid-line of the ring perpendicular to the line of sight. 

The remarkable water maser emission in the nucleus 
of the galaxy NGC 4258 traces a PV diagram where the 
detected emission around the systemic velocity of the 
galaxy comes from the inner edge of the amplifying ring; 
this emission delineates a straight line just as the straight 
line that connects points C and D in Figure 2 (see Figure 
3 of iMivoshi et al.1ll995f ). According to our model re- 
sults, this implies that the absorption coefficient within 
the molecular ring of NGC 4258 is not uniform, instead 
it must be a decreasing function of the distance from the 
central mass, i.e., k = kq^ with q > 1. Moreover, the 
observed red/blue-shifted emission at high velocities that 
arises from the mid-line of the ring perpendicular to the 
line of sight traces a Keplerian curve such as is indicated 
by the model results (see Figure 7). Simply stated, when 
g > 1 the PV diagram is qualitatively similar to the one 
observed for the water maser emission detected in the 
nucleus of NGC 4258. 

As an example, in Figure 7 we show a comparison be- 
tween the results of Model III {q = 15/8) and the water 
maser emission in NGC 4258. The detected emission 
arises from the inner edge of the amplifying ring and the 
mid-line perpendicular to the line of sight. The locus of 
the observed maser emission coincides with the locus of 
the most intense maser emission as indicated by the sizes 
of the circles in Figure 7. The model results indicate that 
there is emission coming from the outer edge of the ring 
at velocities close to 1, nevertheless the sizes of the circles 
indicate that this emission is very weak. Maybe maser 
emission is not detected from this locus for this reason. 

Additionally, our model results also indicate that the 
intense maser emission at the inner edge of the ring ex- 
tends neither to velocities very different from the sys- 
temic velocity nor to impact parameters very different 



from zero, as is indicated by the size of the circles in the 
PV diagram shown in Figure 7. Furthermore, accord- 
ing to the size of the circles, the other locus of intense 
maser emission is the mid-line of the ring perpendicular 
to the line of sight, precisely the locus of the red/blue- 
shifted maser emission at high velocities that describes 
Keplerian curves in the PV diagram. 

4. DISCUSSION AND CONCLUSIONS 

In our model, we have assumed that the gas in the re- 
gion inside the masing ring is transparent to the maser 
radiation. This implies that the most intense maser emis- 
sion at low velocities (velocities near the systemic veloc- 
ity) comes mainly from the outer edge of ring (for q < 1) 
or from the inner edge (for q > 1), either the near or 
far side of the ring, as is indicated in Figure 4. Mea- 
surements of positive acceleration of the maser emission 
around the systemic velocity show that this emission cer- 
tainly comes from the near side of the ring at th e prox- 
imity of its inner edge (e.g.. lGreenhill et al. Ill99"5l ). If we 
suppose that the gas inside the masing ring is thermal- 
ized probably due to its higher density then an important 
fraction of the maser emission from the backside of the 
ring would be absorbed and the detected emission would 
come from the front side of the ring at the outer or inner 
edge depending on the value of q. For instance, consider- 
ing absorption and emission from the gas located inside 
the masing ring, the difference in the intensity for a line 
of sight that passes through both the inner absorbing 
region and the front side of the masing ring from the in- 
tensity for a line of sight that passes through both the 
backside of the masing ring and the inner absorbing re- 
gion is S{1 — e~'^^){e'^^ — 1) where S is the source function 
of the gas inside the masing ring, T2 is the optical depth 
in this region, and ri is the optical depth for the front 
side of the ring. If T2 >> ti, then the detected emission 
would be the radiation amplified by the front side of the 
masing ring. 

Also, we have made a simplifying assumption about 
the geometry of the masing ring in NGC 4258, consider- 
ing that the amplifying ring is strictly flat. Despite the 
observations indicate an apparent warp in th e maser dis- 
tribution of this galaxv. lKartie et al. I ()1999D presented a 
model in which the disk does not require to be physically 
warped in order to the masing gas become exposed to 
the central continuum radiation. In this scenario, dusty 
clouds provide the shielding of the high-energy contin- 
uum, which is required for the gas to remain molecu- 
lar. They found that a flat-disk model of the irradiated 
ring could be applied to a source like NGC 4258 only 
if the water abundance is higher than the value implied 
by equilibrium photoionization-driven chemistry. A very 
important result from their study (based on radiative 
and kinematic considerations) was that, even if the disk 
in NGC 4258 is warped, the maser-emitting gas must be 
clumpy, instead of homogeneous as in the scen ario pre- 
viously proposed by iNeufeld fc Malonev I (|1995f ). 

An important result of our model shows that the as- 
sumption, commonly used, that considers a uniform or 
constant absorption coefficient within the masing ring 
in Keplerian rotation around the nucleus of NGC 4258 
is not appropriate. For example, 'Walli n et al. I (|1998l ) 
supposed that k was constant considering that the locus 
of the maser emission from NGC 4258 was determined 
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mainly by the velocity gradients in a Keplerian velocity 
field indicating some uniformity of k, at least on length 
scales comparable to the coherence or correlation length 
resulting from the Keplerian velocity gradients. On the 
contrary, from our analysis, we conclude that a constant 
absorption coefficient would result in a PV diagram qual- 
itatively different from the observed one, since the most 
intense maser emission would come predominantly from 
a narrow region close to the outer edge of the ring instead 
of a narrow region close to the inner edge of the ring, as 
indicated by the observations. Necessarily, the absorp- 
tion coefficient must be a decreasing function of distance 
from the central mass (i.e., k = kqv^'^ with q > 1) to 
have significant emission coming from the inner edge of 
the amplifying ring and hence explain the form of the PV 
diagram delineated by the water masers in NGC 4258. 

When comparing our edge-on disk model with the 
observations of NGC 4258, it is clear that we need a 
K oc radial dependence for the absorption coefficient 
with g > 1 (so as to favour the emission from the in- 
ner edge of the disk, see above) in order to reproduce 



the observations. In reality, the fact that the disk of 
NGC 4258 is warped introduces geometrical effects which 
might favour the inner disk edge emission (over the one of 
the outer edge). One will need to compute more complex, 
3D transfer models to see whether or not these geomet- 
rical effects are sufficient to explain the PV diagrams of 
the NGC 4258 masers without introducing the radially 
dependent absorption coefficient which is required by the 
edge-on disk models described in the present paper. 
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APPENDIX 

ANALYTICAL APPROXIMATIONS FOR THE OPTICAL DEPTH 

In this appendix, we describe how to obtain the analytical approximations for the optical depth given by equations 
(6)-(8). 

First, we define w = {v — Vr)/^VD then we can change the variables in equation (2) and rewrite it as 



aw 



(Al) 



where 



and 



"^max — 



y-Vr 



Note that Wmm > Wmax since ro < 1. Using the expressions for the line-of-sight velocity v = y/ix^ + y"^)^/^ and the 

r n 1/2 

previously defined variable w = {v — Vr)/^VD, we can write x = {y/{vr + w^vd))^^^ — y^ 

w = we obtain 

2At;^(2/K)2(2-9)/3 



dw ^Vr^/{y/vrY/^ - y2 ■ 



. At zero order around 
(A2) 



Additionally 



exp(— w^)dw = 



erf(wmax) - erf(M;„iin) 



(A3) 



where erf(w) is the error function, defined as erf(ti;) = (2/y^) exp{—t'^)dt. Finally, substituting equations (A2) 
and (A3) into (Al), we obtain the approximation for the optical depth 

20FKoAi;z5(y/«,)2(2-9)/3 



^{Vr,y) 



erf(wmin) - erf(Wmax) 



(A4) 



Around the inner edge of the ring, Vr — yl'i'^J'^ and the maximum value of [erf(u>min) — crf(?i>n 



= 2. Substituting 

these approximations into equation (A4), we obtain equation (6). Similarly, around the outer edge of the ring, Vr — y, 
and the maximum value of [erf(winin) — erf(wmax)] also equals 2. Then, equation (A4) reduces to equation (7) for the 
optical depth at the outer edge of the ring. 

In order to find an approximation for the local maximum of optical depth at the mid-line of the ring perpendicular 
to the line of sight, we expand the expression for the velocity along the line of sight around a; = to obtain 

1 3 x^ 
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therefore 



and thus 



V — 



,1/2 



3 



= -Au 



D ■ 



^ ~ 31/2 ^^£> i' 



Then, using Sobolev's approximation 



r = Ko{x^+yY''^\2x), 
and substituting equation (A7) into (A8), we obtain the following expression 

since y <C 1 and Avd is small, then ^AvdV^^^ ^ 2/^) hence 

3I/2 i3 » ' 

considering that y — v^'^, we finally obtain the approximation for the local maximum of the optical depth at the 
mid-line of the ring given by equation (8). 
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Fig. 1. — Schematic diagram of a gaseous disk in Keplerian rotation. Tlie masiiig gas exists between radii Rq and R. At radii smaller 
than Ro and greater than R. the disk is transparent to the maser radiation. Tlie observer is on the plane of the disk. All the distances axe 
measured in units of R, the outer radius of the amplifying ring; therefore the variables x, y, and ro are dimensionless. 
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Fig. 3. — Prom top to bottom results of Models I, II, III, and IV. The left panels show PV diagrams for the maser emission peak. The 
filled squares, circles, and triangles represent the strongest maser emission which is coming from regions o, b, and c, respectively, indicated 
in Figure 4. The straight lines or curves represent the velocity dependences of the regions where this emission arises. The central panels 
show PV diagrams for the maser emission peak. Because of the point-symmetric shape of these diagrams, only positive velocities are shown. 
The radii of the open circles are proportional to the maximum maser intensity at each position and velocity. The right panels show the 
logarithm of the ratio between the maximum intensity and the background intensity as a function of the velocity. 
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Models I and II Model IE Model IV 




Fig. 4. — Schematic representation of the results of Models I, II, III, and IV. The filled squares, circles, and triangles represent the most 
intense maser emission that is coming from regions a, b, and c, respectively. These regions are very narrow because the correlation distance 
is very small; that is, the width Avd is much smaller than the line-of-sight velocity gradient. Besides, the exponential amplification of the 
intensity emphasizes small changes in the optical depth. 




Fig. 5. — Left: PV diagrams for the mascr emission in grey scale with intensity contours overlaid. The darker regions show the locus 
of the strongest emission in these diagrams, which potentially could be detected, depending on the sensitivity cutoff of the observations. 
Right: Close-up to the PV diagrams showing the maser emission at low velocities. 
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Fig. 6. — Vc as function of the inner radius of the ring, ro- Bottom: For Vc < 1, it is computed from equation (9a) for some representative 
values of q. For g < 1, there is no solution to equation (9a). Top: For Vc > 1, it is computed from equation (9b). This is the solution to 
equation (9b), using Avu = 0.007t)out and q = 5. 
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Fig. 7. — Comparison between the calculated PV diagram for the maser emission peak in Model III {q = 15/8) and the PV diagram 
deUneated by the water masers observed in NGC 4258 (Miyoshi et al. 1995). The radii of the open circles are proportional to the maser 
intensity at each position and velocity. The dots represent the observed maser spots in NGC 4258. We have subtracted the ring systemic 
velocity of 476 km from the observed local standard of rest velocity of the maser spots in order to compare the observed PV diagram 
with the modeled one. The positions and velocities are in units of the outer radius of the ring (8 mas) and the rotation velocity at the 
outer edge of the ring (770 km s~^), respectively. 



